Eigenfunction fractality and pseudogap state near superconductor-insulator transition 
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We develop a theory of a pseudogap state appearing near the superconductor-insulator transition 
in strongly disordered metals with attractive interaction. We show that such an interaction combined 
with the fractal nature of the single particle wave functions near the mobility edge leads to an 
anomalously large single particle gap in the superconducting state near SI transition that persists 
and even increases in the insulating state long after the superconductivity is destroyed. We give 
analytic expressions for the value of the pseudogap in terms of the inverse participation ratio of the 
corresponding localization problem. 



Rapidly growing number of experimentsii 2 ^ 3 ^^^^ 
on various disordered superconductors show that a 
novel phase often appears on the insulating side of the 
superconductor-insulator transition. As the disorder 
strength is increased the superconductivity is suppressed 
leading to a strange insulator characterized by a large 
thermally assisted resistance with a small but hard gap 
(Fig. 1). Further increasing the disorder one gets the 
usual variable range hopping behavior. Experimentally, 
the phase diagram of disordered superconductors is often 
explored by varying the applied magnetic field. On the 
superconducting side of the transition a relatively small 
field destroys the superconductivity resulting in a hard 
gap insulating state. At larger magnetic fields the resis- 
tance and the gap drop^SiL^. This is observed only 
in a narrow window of disorder strengths, away from 
this window on a superconducting side the application 
of magnetic field converts the supercondcutor into a nor- 
mal metal as usual. 

It is tempting to explain these data by the formation of 
localized Cooper pair o 4 i 10 . In this picture the supercon- 
ductivity is due to a fragile coherence between localized 
Cooper pairs, while the energy to break the pair is much 
larger and remains finite even when the coherence (and 
thus the superconductivity) is destroyed. The hypoth- 
esis of preformed Cooper pairs is further confirmed by 
the behavior of these superconductors at higher temper- 
atures. On the insulating side of the transition in thick 
(effectively 3D) films one observes 1 ^ Arrenihus temper- 
ature behavior of the resistivity, R(T) cx exp(T//T), at 
low temperatures. The experimental value of the activa- 
tion energy, Tj, is somewhat larger than the supercon- 
ducting gap in less disordered samples and grows with 
the disorderii However, at higher temperatures this 
behavior is replaced^ by Mott's variable range hopping 
R(T) ~ exp(T M /r) 1 / 4 . This can be understood if the 
insulating pseudogap is due to preformed Cooper pairs 
with a relatively large pairing energy Tj. 

In this Letter we show that preformed Cooper pairs 
appear in the semi-microscopic three-dimensional model 
that contains only low energy electrons with weak BCS- 
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FIG. 1: Schematics of low-temperature phase diagram 
of disordered superconductors in the vicinity of the 
superconductor-insulator transition as a function of the mag- 
netic field, B, and disorder strength. 



type attraction and a strong random potential that leads 
to Anderson localization of single-particle states. 

In the presence of preformed Cooper pairs the parity 
effect should arise - the ground state energy with even 
number of electrons is lower then that for the closest 
odd number. The corresponding parity gap Ap ~ Tj ~ 
5/ ln((5/A) has been calculated by Matveev and Larkin 11 
for small superconducting grains. Here 5 — l/(v L 3 ) is 
the mean level spacing in the grain, A << S is the energy 
gap in the bulk superconductor, and vq and L 3 are the 
density of orbital states and the volume of the grain, 
respectively. In this Letter we argue that the result of 
Refiii can be generalized to bulk Anderson insulators. In 
this case, L is replaced by the localization radius L\ oc , so 
that 5 — » 5l = l/(^o£ 3 oc ), and the BSC A is replaced by 
the superconducting gap A cr ;t at the Anderson transition 
point. The fractal nature of near-critical wavefunctions 
characterized by the fractal dimension D2 < 3 determines 
the reduction factor (<5i/A cr it) 1_D2 ^ 3 S> 1 that replaces 
ln((S/A) in the Matveev-Larkin formula. 
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We consider two different regimes. In more disordered 
materials where 5l 3> A cr ; t the Cooper instability and 
superconductive long-range order disappear. However, 
the attraction between electrons persists as long as 5p 
remains smaller than the Debye frequency ujb and results 
in the "local" pairing of electrons with opposite spins 
occupying the same localized state. We show below that 
in this regime the hard-gap insulator is formed whith 
properties similar to those observed ini£. We further 
argue that in the less localized regime, when bp < A cr ; t , 
the unusual superconductive state with a pseudogap is 
formed. The features of this state are: (i) single-electron 
excitation gap Ai is larger than superconductive gap A 
so that the ratio Ai/T c is anomalously high, and (ii) 
insulating trend in the R{T) curves exists above T c . 

We assume that superconductivity is due to attraction 
between electrons that originates at high energy scales 
~ u>d and that it is affected by localization of electron 
wave functions only for a very large disorder Sp > LOp>. In 
a fcrmion system with weak attraction one can leave only 
the pair interaction terms in the Hamiltonian leading to 
the usual BCS model in the basis of localized electron 
states 12 : 



H = 
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where 



M 



drtf(r)rt(r), 



(1) 



A is dimensionless Cooper coupling constant, tj is the 
single-particle energy of the state j, and cy CT is the corre- 
sponding electron operator for the spin projection a. 

Physical properties of the electron system are con- 
trolled by the electrons near the Fermi level, so a very 
important implicit ingredient of the model Q is the 
statistics of matrix elements Mjk between eigenstates in 
the vicinity of Anderson mobility edge. The key fea- 
ture of these nearly-critical wavefunctions is their frac- 
tal structure^ that shows in the anomalous scaling of 
diagonal matrix elements Mjj = Mj ("inverse partici- 
pation ratios", IPRs) with localization length: typical 
IPR M cx i;~f 2 , where fractal dimension D 2 < 3. Nu- 
merical studiesi^i^ indicate that Di = 1.30 ± 0.05 for 
the standard 3D Anderson transition. The IPR distribu- 
tion function V(Mj) has been studied ini^. Scaling the- 
ory of localization predicts that near the mobility edge 
V(Mj) acquires a scale- invariant form and this is indeed 
what was observed^. The same data demonstrate that 
V(Mj) decreases fast for atypicallly extended states, i.e. 
at Mj/M <c 1. This allows us to use the average value 



M = Lq 3 (L 1oc /L q ) 



-D 2 



(2) 



where Lq is the short-scale cutoff length of the fractal 
behaviour. The associated energy scale En = 1/{vqLq) 
depends on the microscopic details of the model of dis- 
order and might be small compared to Fermi-energy Ep. 
Localization length depends on Fermi-energy (in the scal- 
ing region L ioc > L ) as L[ oc w L (E /(E m - E F )) V , 



where E m is the position of the mobility edge and v is 
the localization length exponent. 

Another important property of the nearly-critical 
eigenstates is their strong correlation in energy and real 
space^£^£^£ even in the limit of strong fractality D2 << 
3. It results in the scaling dependence of the average ma- 
trix elements Mjk on the energy difference Ej — E^ = lj: 



VM 



M(w) 



(Li oc /L a ) 3 t at w<(5 L (a) 
(E /cj)-y at S L < uj < E (b) 



(3) 

where 7 = 1 — D2/3 and V is the total system's volume. 
Note that in the critical region M(u>) ^> 1 in contrast 
both to a metal and to a deep insulator (Li oc ~ La), 
where M(gj) w 1. 

We begin with the insulating region A 5l *C u>n 
where Cooper interaction can be treated perturbatively. 
In the first order of the perturbation theory we take into 
account only diagonal terms j = k of the interaction 
similar to the case of ultrasmall grainAi. Then the energy 
(counted from Ep) required to break a bound pair of 



ii) 



electrons siting in the j-th orbital state is 2Ap 
Typical value of this " parity gap" (cf>ii) scales then as 



A, 



-En 



L 
L\ oc 



D-2 



cx (E m - Ep) 



vD 2 



(4) 



Neglecting fluctuations in the local values of Ap', one 
finds that all states occupied by single electrons are 
shifted up by the amount Ap that leads to the electron 
DoS v(e) = v a d{e - A P ). In fact, values of Ap* differ 
for different localized states, and the average density of 
states v(e) in a large sample is determined by the IPR 
distribution V{M): 



v{e) = v n 



2ev /\ 



V{M)dM 



(5) 



As mentioned above, numerical data for V(M) indicate 
its very fast decrease at M/M — * 0. Thus the DoS 
shape (JSJ is not far from a rectangular sharp gap, with 
the gap value given by Eq. We emphasize that i) 
parity gap Ap is much larger than level spacing 8l at 
Lioc/Ln > (2/A) 1 ^ 1 ^", and ii) the DoS © does not con- 
tain any " coherence peak" above the gap (cf. Ref. ITij) . 

We associate^ the spectral gap Ap with the mea- 
suredi^ activation energy Tj. The external parameter 
(E m — Ep) representing the disorder strength in Eq. Q 
can be replaced with an experimentally more accessible 
parameter (er c — a) cx (E m — Ep). Here a is the high 
temperature conductivity and o~ c is the value of the con- 
ductivity where the parity gap Ap first develops. We 
obtain 



Ti = A(l - a/o- c y l 



(6) 



where A is conductivity-independent. This equation pre- 
dicts a moderate increase of Tj with disorder strength in 
agreement with the experimental datai, see Fig. 2. 
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FIG. 2: Experimental values of the gap from RefQ, Tj (boxes) 
and a fit to the equation © with v = 1, D2 = 1.3. The only 
fitting parameter was the constant A = (A/2)i5o; the data 
points of Ref. [1] correspond to Eo ~ WOK at A ~ 0.2 ex- 
tracted from the BCS value of T c ~ 3K for less disordered 
samples- and ujd ~ 500K. The value of a c was determined 
from high T data. Application of scaling formulas to the 
data shown here is justified by the large value of the local- 
ization length of the most disordered sample shown in this 
plot: L^ x > 30A which was deduced from the Mott tem- 
perature characterizing the resistivity of similar samples at 
intermediate temperaturesi- 



We now turn to the parameter region < A where 
one expects a global superconductive coherence to exist 
at low enough T. Indeed, in this regime a given localized 
single-particle state typically overlaps in a real space with 
a large number ~ A/Sl of eigenstates tpj in the same 
energy strip |ej| < A. It is natural to expect that in 
this case the mean-field approximation can provide some 
guidance, while not necessarily being quantitatively ac- 
curate (e.g. because fluctuations of the overlap matrix 
elements Mjk are large). To test the validity of the mean- 
field, we compared its prediction for T c (Sl) (see below) 
with the transition temperature that was found numeri- 
cally by computing the first terms of the virial expansion 
applied to Q with My determined by exact diagonal- 
ization. Reasonably good agreement was found-S. 

To proceed with the mean-field analysis, we introduce 
averaged energy-dependent pairing amplitudes F(ej) — 
(cj^Cji) and the gap function A(e) = A f deiM(e — 
ei)F(ei). Following standard steps, we decouple inter- 
action term in the Hamiltonian Q via the gap function 
A(e), calculate anomalous averages F(e), and arrive at 
the modified BSC gap equation in the form 



A(e) 



xAi + A 2 ( £l ) 



2T 



(7) 



Gap function A(e) obeying Eq.(J7| is an even function of 
e with the maximum value A(e = 0) = Ao. 

Superconducting transition temperature T c is deter- 
mined by linearization of Eq.|7|) with respect to A(e). 
Due to power-law decrease of M(tu) at large arguments, 
the integral in Eq.Q converges and is dominated by 
lj ~ T ("infrared superconductivity") so no upper cut- 
off is needed, contrary to usual BSC problem. When 



Fermi level is very close to the Anderson mobility edge 
E m and level spacing Sl ~ E ((E m - E F )/E a ) 3ly is neg- 
ligibly small, one can use for M(u>) the Eq.J^)). Then 
the critical temperature is given by 

r c °(A, 7 )= J BoA 1 ^C( 7 ) (8) 



where dimensionless function C{^) can be computed nu- 
merically. At small A the value given by Eq.© exceeds 
the BCS value Tbsc ~ ud^ 1 ^- This may lead to a 
maximum in T c near the critical disorder. The zero- 
temperature energy gap At=o m the same limit Sl — > 
is given by Eq.© with C(j) replaced by another function 
D{n). We plot D^) and 2£>(7)/C(7) in Figure 3. 
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FIG. 3: Suppression of T c as a function of level spacing within 
localization volume. The insert shows 7 -dependence of the 
dimensionless gap and 2A(0)/T C for Sl = 0. 

Using Eq.© we eliminate the interaction constant A 
and the cut-off parameter Eq from Eq.© and arrive at 



2(C7( 7 ))i 



(T°Jh 



(9) 



This formula (applied at Sl 3> T®) generalizes the re- 
sults of Ref^i to bulk strongly disordered superconduc- 
tors. Note that in contrast to Refi^i here the reduction 
of Ap compared to Sl is not due to the renormalization 
of attractive interaction-i, but to the enhancement of the 
matrix elements Mj due to fractality. 

To study the effect of nonzero level spacing Sl upon T c , 
we approximate M{uj) by a simple interpolation formula 
M{ui) = E^ito 2 + 5 L )~i and solve the linearized version 
of Eq.© for T c {5l) numerically. Since M(w) is a uni- 
form function of T and Sl , while the coupling constant 
A and Eq enter Eq. © only in a combination EqX 1 ^ , it 
is possible to present the dependence T c (<5l) in the form 



T..iS L ) = T!!T, ('J 



(10) 



where Xj? is defined in Eq.© and scaling function 7^,(x) 
does not depend on A. This universal function was found 
numerically for 7 = 0.57 (corresponding to D2 = 1.3 for 
3D Anderson transition), the result is plotted in Fig. 3. 
Function Xy(x) approaches very slowly its asymptotic 
value 1 at large x. The actual T c is suppressed as com- 
pared to the mean-field result due to fluctuations. We 
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expect it still has the scaling form i|10|l but with a differ- 
ent scaling function that vanishes at x < x c ~ 1. 

Although gap equation Q is similar to the conven- 
tional one, the real-space properties of the state that 
it describes are unusual. The local pairing amplitude 
F{ Y ) = Si'( c iT c jl)'^'i ( r ) i s extremely inhomogeneous in 
space, populating only small fraction cx [T c (5 l) I EqV of 
the total volumeSS. Diamagnetic response of such a su- 
perconductor differs strongly from that of usual "dirty- 
limit" materials with uniform |_F(r)|, but reminds that 
of weakly-coupled Josephson junctions arrays. Qualita- 
tively, we expect (i) extremely weak Meissner effect and 
considerably stronger linear shielding effect, and (ii) su- 
perconductor — > gaped insulator — > Mott insulator se- 
quence of transitions upon magnetic field increase. 

We now discuss the effect of "local pairing" consid- 
ered previously (see Eq.@J) in insulating state. We have 
seen that in insulating state when 5l 3> A single-particle 
excitation (carrying spin i) have a gap Ai = Ap. On 
the other hand, excitations that involve only hopping of 
paired electrons between localized levels and do not in- 
volve breaking pairs are gapless, i.e. the energy gap A 
vanishes. On the other hand, in superconducting state 
within mean-field we have Ai = A + Ap, while the gap 
for pair excitations^ (without pair breaking) is A2 = 2A. 
When Li oc /Lq — > 00, the parity gap Ap becomes much 
larger than Sl- Therefore, we expect that there is a 
regime Ap > A in the superconducting state where the 
spin gap Ai is larger than the energy gap. 

In the derivation of the gap equation (JJJ) we neglected 
the effects of the parity gap which is correct for Ap <C T. 
A large parity gap pushes up all unpaired states with 
respect to the paired ones which modifies the thermal 
distribution factor in the Eq.0. In the limit of Ap 3> T 
the contribution of unpaired electrons vanishes and the 
problem becomes equivalent to the pseudospin models. 
This effectively doubles the value of the gap and thus 
replaces 2T — + T in the tanh() argument in Q. For a 
general Ap ~ T the crossover formula between these two 
regimes can be obtained in the mean-field approximation. 
The result is (for derivation cfi^fi ) that tanh(...) in Eq.Q 
should be replaced by 



sinh B 



cosh B 



,-Ap/T 



where B = V(e? + A 2 (ei)/T, and e x = e x - A P . 

Note that most experiments probing the superconduct- 
ing gap, such as tunnelling conductance, optical conduc- 
tivity or NMR, measure the spin gap Ai. One thus 
expects that these data would show anomalously large 
(compared to BCS value of 1.76) ratio of spectral gap to 
transition temperature near S-I transition. An additional 
suppression of T c in comparison with Ap can be due to 
electron-electron interaction in the density channel (ir- 
respectively of the interaction sign) that is not included 
into the model |QJ. 

An anomalously large ratio A\jT c leads to the insulat- 
ing trend of the resistivity versus temperature behavior in 



the intermediate temperature range T c < T < A x . This 
was observed in strongly disordered superconductors, see 
e.g£, and is well-known as a pseudogap phenomenon in 
underdoped cuprates^*. The quantitative similarity be- 
tween R(T,B) behavior in InO^ films and underdoped 
cuprate o 2 ?i 2 ? allows one to speculate that the pseudogap 
in underdoped cuprates might also be related to pairing 
of electrons on localized states. The important difference 
of the cuprates is the d-wave symmetry of the pairing. 

The effect of magnetic field on the fractal superconduc- 
tor deserves a separate study. Here we only remark that 
the available magnetoresistance data are not consistent 
with a conventional picture of weakly coupled supercon- 
ducting grains but point towards fractal superconductiv- 
ity. Typically) 4 i^i7i?, the magnetoresistance is negative 
at high fields and persists up to very high field values. 
Consider, for instance, the magnetoresistance data on 
amorphous InO^, shown in Fig.l of Ref. U (the samples 
were of the same type and origin as those studied ini*£) . 
Activation energy for this sample is Tj ss 15 K (as deter- 
mined in the temperature range 1.3 — 5 Kelvin); applica- 
tion of 8 Tesla field at the temperature of 1.58 K leads to 
negative magnetoresistance (MR) R(H = 8T)/i?(0) « 
0.5. Interpreting this effect as being due to magnetic- 
field suppression of the pairing gap in small grains, we 
find Ti - Ti(H = 8T) « 0.7 K, which is about 5% of Tj 
only. Interpolating this dependence we estimate that the 
field necessary to destroy completely superconductivity 
in each grain is about H^ p ~ 50 — 80Tesla™ which is 
much larger than one expects for realistic grain sizes. In- 
deed, orbital critical magnetic field for a small (radius 
R < £, where £ = WHD/A is the coherence length) su- 



26,27 



to H%f 



1 000 T 



where 



perconductive grain is equal 
R and £ are measured in nanometers. Using a typical 
diffusion constant for a poor metal D « lcm 2 /s and 
allowing for a very high gap value A = 10 K, we find 
£ = 8.5 nm. Together with the lowest bound for the 
grain radius R = 6 nm (determined via the condition 
5 = (A^R 3 )' 1 < A) it leads to H e c f < 20 Tesla which 
is still much smaller than H°* p above. These estimates 
did not take into account the spin effect of magnetic field 
that would further decrease the value of . This dis- 
crepancy between the experimental value of the char- 



(11) acteristic magnetic field, ff^ p , and its estimate, 



can be accounted for by the fractal nature of the paired 
states because the effect of magnetic field on the local 
pairing gets much weaker due to low value of fractal di- 
mension Z?2 = 1-3 of these states that indicates "almost 
one-dimensional" nature of these eigenstates. 

In conclusion, weak Anderson insulators with Cooper 
attraction are shown to possess hard insulating gap 
whose magnitude is determined by the IPR statistics 
near the mobility edge. Although this gap is due to 
electron pairing, it does not lead to a coherence peak. 
In the ground-state of this insulator all electrons are 
paired on individual localized eigenfunctions. When the 
Fermi-level gets closer to the mobility edge, superconduc- 
tive correlations develop between localized pairs. Key 
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features of the predicted superconductive ground-state 
are extreme inhomogenity of superconductive correla- 
tions in real space, an unusually large (compared to T c ) 
single-particle excitation gap (spin gap), and pseudo- 
gaped regime at temperatures about T c . All these un- 
usual features are due to the fractal nature of localized 
eigenstates near the mobility edge. 
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